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ABSTRACT. A one-dimensional Hamiltonian system with exponential interactions per- 
turbed by a conservative noise is considered. It is proved that energy superdiffuses and 
upper and lower bounds describing this anomalous diffusion are obtained. 

1. Introduction 

Over the last decade, transport properties of one-dimensional Hamiltonian systems con- 
sisting of coupled oscillators on a lattice have been the subject of many theoretical and 
numerical studies, see the review papers [9, 1 1, 23]. Despite many efforts, our knowledge 
of the fundamental mechanisms necessary and/or sufficient to have a normal diffusion re- 
mains very limited. Nevertheless, it has been recognized that conservation of momentum 
plays a major role and numerical simulations provide a strong evidence of the fact that one 
dimensional chains of anharmonic oscillators conserving momentum are superdiffusive. 

An interesting area of current research consists in studying this problem for hybrid mod- 
els where a stochastic perturbation is superposed to the deterministic evolution. Even if the 
problem is considerably simplified, several open challenging questions can be addressed 
for these systems. In [2] it is proved that the thermal conductivity of an unpinned harmonic 
chain of oscillators perturbed by an energy-momentum conservative noise is infinite while 
if a pinning potential (destroying momentum conservation) is added it is finite. In the same 
paper, diverging upper bounds are provided when some nonlinearities are added. This does 
not, however, exclude the possibility of having a finite conductivity. Therefore much more 
interesting would be to obtain lower bounds showing that the conductivity is infinite and 
that energy superdiffuses, but this problem is left open in [2]. 

In [7], has been introduced and studied numerically, a class of Hamiltonian models for 
which anomalous diffusion is observed. There, the investigated systems present strong 
analogies with standard chains of oscillators. They can be described as follows. Let V and 
U be two non-negative potentials on M and consider the Hamiltonian system ( r (t ) , p(f ) ) f >o 
whose equations of motion are given by 

^=V(r x+1 )-V(r x ), ^=U'(p x )-U'(p x+1 ), x G Z, (1) 

where p x is the momentum of the particle x, q x its position and r x = q x — q x _ \ is the 
"deformation" of the lattice at x. Standard chains of oscillators are recovered for a quadratic 
kinetic energy U(p) = p 2 /2. Now, take V = U, and call Ij^c-i = f x and rj2 x = p%- The 
dynamics can be rewritten as: 

dru(f) = (y'(r] x+l )-V'(n x ^)yt. (2) 
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Notice that with these new variables the energy of the system is simply given by Y.xeZ^i'Hx)- 
In [7] an anomalous diffusion of energy is numerically observed for a generic potential V. 
Then, following the spirit of [2], the deterministic evolution is perturbed by adding a noise 
which consists to exchange rj x with 77 T+ i at random exponential times, independently for 
each bond {x,x + 1}. The dynamics still conserves the energy V(J7. T ) and the "vol- 
ume" T,xez r l^ an d destroys all other conserved quantities. As argued in [7], the volume 
conservation law is responsible for the anomalous energy diffusion observed for this class 
of energy-volume conserving dynamics. This can be shown for quadratic interactions ([7]) 
with a behavior similar to the one observed in [2]. For nonlinear interactions the problem 
is much more difficult. 

The aim of this paper is to show that if the interacting potential is of exponential type 
then the energy superdiffuses. Therefore, for this class of related models, in a particular 
case, we answer to the open question stated in [2] . With some additional technical work 
we think that our methods could be carried out to the Toda lattice perturbed by an energy- 
momentum conserving noise (considered e.g. in [19]). 

The paper is organized as follows. In Section 2 we define precisely the model. The 
results are stated in Section 3. To prove the theorems we first perform a microscopic 
change of variables (Section 4) which permits to use a nice orthogonal decomposition of 
the generator (Section 5). Roughly speaking the upper bound on the energy superdiffusion 
is proved in Section 6 and the lower bound in Section 7. Section 8 contains a comment 
about the possible extensions and comparisons of our model to others. In the Appendix we 
prove the existence of the infinite dynamics. 

Notations: For any a,b £ M 2 , a ■ b stands for the standard scalar product between a 
and b and \a\ = s/a ■ a for the norm of a. The transpose of a matrix M is denoted by M T . 
If u : x = [x\ , . . . ,x„) € W — s- m(x) = (u\ (x), . . . ,m^(x)) e Mr is a differentiable function 
then (d Xj Uj)(x) denotes the partial derivative of uj with respect to the i-th coordinate at x 
and (Vm)(x) denotes the differential matrix (the gradient if d = 1) of u at x, i.e. the d x n 
matrix whose (i,j)-th entry is (d Xj Uj)(x); if u := (mi, . . . ,u^) : Z — > W 1 then we adopt the 
same notation to denote the discrete gradient of u defined by (Vm) := (Vwi, . . . , Vm^) with 
(Vui)(x) = iii(x+ 1) — ui(x). 



Let b > and Vb(q) = e~ bq — l+bq. We consider the system 77(f) = {T] x (t) : x e Z} 
on R z defined by its generator L = A + yS, J > 0, where for local 1 differentiable functions 
/ : R z -> R we have that 



The deterministic system (2) with potential Vb is well known in the integrable systems 
literature. It has been introduced in [20] by Kac and van Moerbecke and was shown to be 
completely integrable. Consequently, the energy transport is ballistic ([7, 30]). As we will 



A function / defined on an infinite product space is said to be local if it depends only on its variable through 
a finite number of coordinates. 
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see this is different when the noise is added: the energy transport is no more ballistic but 
superdiffusive. 

The existence of the dynamics generated by L is proved in the Appendix for a large set 
of initial conditions and in particular for a set of full measure w.r.t. any invariant state fip ^ 
(see bellow for its definition). 

The system conserves the energy YaeZ^bi'Hx) and the volume £xeZ 'Hx- In f act > we have 

where the microscopic currents are given by 

7x,x+i (J?) = -b 2 e- b ^ + ^ +b 2 (e- b ^+e- h ^)- /V(V^)) 

and 

Lx+x (n) = be- b ^+be- b ^ - 7V77,. 
Every product probability measure ^ on M. z in the form 

HA ^ ) = II 2 1 (P ~ X ) ex P{ -Pe~ bTk - A J?,} , P > , A > 

xeZ 

is invariant for the dynamics. 

Let denote the average with respect to ju^ We define e := e(j3,A),v := v(j3, A) 
as the averages of the conserved quantities Vb(rj x ), rj x with respect to /ig ^, respectively, 
namely e = (V b (rj x )} ^ and v = (j?x) M(U ■ 

A simple computation shows that 

(Jx*+i)n ?iX =b 2 (e-bv) 2 and (]' xx+x ) ^ 1 =2b{e -bv). (3) 

Hence, in the hyperbolic scaling, the hydrodynamical equations are given by 

d t t-b 2 d q ((t-b») 2 )=0 
d t t) + 2bd q (e-bx>) = 

and can be written in the compact form d t X + dq3(£) = with 

*-(:)• )■ 

This can be proved before the appearance of the shocks (see [7]). The differential matrix 
of 3 is given by 

M)(S)-2» (-*-*> "C-")). 

For given (c,v) we denote by (f r + ) r >o (resp. (r r ~) r >o) the semigroup on S(R) x 5(R) 
generated by the linearized system 

d,e+M T d q e = 0, (resp.d t e-M T d q e = 0), (6) 

where 

M:=M(e,v) = [V5](ffl), <» = 

We omit the dependence of these semigroups on (e, v) for lightness of the notations. Above 
S(M) denotes the Schwartz space of smooth rapidly decreasing functions. 




4 



CEDRIC BERNARDIN AND PATRICIA GONCALVES 



3. Statement of the results 

2 d z _ji 

For each integer z > 0, let H z (x) = (— Vfe be trie Hermite polynomial and 

h z (x) = (zl\/27z)~ l H z (x)e~ x2 the Hermite function. The set {h z ,z > 0} is an orthonormal 
basis of L 2 (K). Consider in L 2 (]R) the operator Kq = x 2 — A, A being the Laplacian on R. 
For an integer k > 0, denote by Hj the Hilbert space induced by 5(K) and the scalar product 
(•, defined by (f,g)k = (/>^og)o> where (•, -)o denotes the inner product of L 2 (M) and 
denote by H_j the dual of H^, relatively to this inner product. Let (■) represent the average 
with respect to the Lebesgue measure. 

We take the infinite system at equilibrium under the Gibbs measure jJ.^ ^ corresponding 
to a mean energy e and a mean volume v. Our goal is to study the energy-volume fluctuation 
field in the time-scale tn l+a , a > 0: 

9f> a (G) = 4= £ G(x/n) ■ (cD,(f« 1+a ) - Co) , (7) 



where for ^el.ieZ, 



<?2(?) ; ' V 

and Gj , G2 are test functions belonging to S(M). 

If £ is a Polish space then D(M + ,£') (resp, C(R + ,.E)) denotes the space of Zs-valued 
functions, right continuous with left limits (resp. continuous), endowed with the Skorohod 
(resp. uniform) topology. Let Q" a be the probability measure on D(R + ,H_,(. x 
induced by the fluctuation field and jip ^ . Let - denote the probability measure on 

Z)(R + ,R Z ) induced by (f](f))r>o an d /^g \- Let E^ - denote the expectation with respect 

Theorem 1. Fix an integer k > 2. Denote by Q the probability measure on C(ffi. + ,]HI_£ x 
H_fc) corresponding to a stationary Gaussian process with mean and covariance given 
by ' 

E e [^(H)^(G)] = (frH.*f-G) 

for every < s < t and H, G in x H*. //ere J := j(j3 , A) ;s f/ze equilibrium covariance 
matrix 2 of (do- Then, the sequence (Q"'°) n >i converges weakly to the probability measure 

Q. 

As a consequence of Theorem 1 we obtain a Central Limit Theorem for the energy flux 
and for the volume flux through a fixed bond. For that purpose fix a site x £ Z, let <? x " x+1 (t ) 
(resp. y" , j (f)) denote the energy (resp. volume) flux through the bond {x,x+ 1} during 
the time interval [0,fn], defined as: 

C+i(0:= I {V b (r]y(tn))-V b (7] y (0))} (resp. *£ +1 (f):= £ {^(/n) -Tj,(0)} 

y>x+l y>x+l 

Then, from the previous result we conclude that 

Corollary 1. Fix x £ Z and let Z"' e := ^{^"x+i (t) - E M ^ [<^ +1 (/)]}. For every & > 
1 and every < t\ < 1 2 ■ ■ ■ <th (Z,"' 6 , • • • ,Z^' e ) converges in law to a Gaussian vector 



See (18) for an explicit value. 
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(Z^ , • • • ,Zf ) with mean zero and covariance given by 

E Q [Z?Z e s ] = ^-(X-bP) 2 s, 

for all s <t. 

Corollary 2. Fix xeZ, let Z r "' v := - % A [^ +1 (f)]}. For every * > 

1 and every < t\ < ti ■ ■ ■ < t^, (Z'"\ ■ ■ ■ ,Z^' 1 ') converges in law to a Gaussian vector 
(Z v t{ , • • • ,Zy ) with mean zero and covariance given by 

e q [z?z v s ] = is, 

for all s <t. 

We notice that, according to Corollary 1, the limiting energy flux Zf has a vanishing 
variance for A = bj5 which is equivalent to e = bv. 

The theorem above means that in the hyperbolic scaling the fluctuations are trivial: 
the initial fluctuations are transported by the linearized system of (4). To see a nontrivial 
behavior we have to study, in the transport frame, the fluctuations at a longer time scale 
tn l+a , with a > 0. Thus, we consider the fluctuation field ( 3/ nM ', a > 0, defined, for any 
G e S(M) x S(R), by 

C"'"(G) = % n a {T+aG) . (8) 
According to the fluctuating hydrodynamics theory ([26], pp. 85-96), in the case of a 
normal (diffusive) behavior a = 1, the field {W"' a ) t >o should converge to the stationary 
field (@t)t>o simply related to the solution (£% )i>o of the linear two dimensional vector 
valued (infinite-dimensional) stochastic partial differential equation 

a3 = v-(^v^) + ^/2^fv-w f . (9) 

Here W(x,t) is a standard two-dimensional vector valued space-time white noise and the 
coefficient £P := 2$(e, v) is expressed by a Green-Kubo formula ( see (12)). 

Our second main theorem shows that the correct scaling exponent a is greater than 1/3: 

Theorem 2. Fix an integer k > 1 and a < 1/3. Denote by Q the probability measure 
on C(R + ,H_£. x H_£.) corresponding to a stationary Gaussian process with mean and 
covariance given by 

E e [^(H)^(G)] = (H -xG) 

for every < s < t and H, G in Hjfc x H^. Then, the sequence (Q"' a ) n > l converges weakly 
to the probability measure Q. 

As in the hyperbolic time scale from the previous result we obtain limiting results for 
the energy and volume flux. In this case, since we take the system in a time dependent 
reference frame, we defined the energy and volume flux through the time dependent bond 
{u*' a (n),u*' a (n) + 1}, where u*' a (n) := —2bX/[}tn l+a . For that purpose, fix a site x € Z 
and let $ x,a, > (resp. Y'l. a , At)) denote the energy (resp. volume) flux through the bond 

{u x t ' a {n),u x t ' a {n) + 1} during the time interval [0,f« 1+a ], defined as: 

#« W W:= E {n(r7,( f « 1+a ))-V fc (r7,(0))}, 

v>«, (n) 

/ \ (10) 

(resp. *£, (B) (f):= £ {rfc(tn 1+a )-Tfc(0)}). 

y>u? a (n) 



(1 
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Then, from the previous result we conclude that 



Corollary 3. Fix t > 0, x e Z and a < 1/3. Then 




0. 



anc/ 



lim E 



0. 



Similar results have been obtained in [18] by one of the authors for the asymmetric 
simple exclusion process. The proof of the three corollaries above follows the same argu- 
ments as in [ 1 8] once the previous theorems are proved. For that reason we will only give a 
sketch of their proof. The proof of the theorems is more problematic since the multi-scale 
analysis performed in [18] relies crucially on the existence of a spectral gap so that we 
cannot follow [18]. Therefore we propose an alternative approach based on computations 
of some resolvent norms. 

Theorem 2 does not exclude the possibility of normal fluctuations, i.e. a = 1. In order 
to show that the system we consider is really superdiffusive we will show that the transport 
coefficient which appears in (9) diverges so that the correct scaling exponent a is strictly 
smaller than 1 . Our third result, stated bellow, shows it is in fact less than 3 /4. 

With the notations introduced in the previous section, the normalized currents are de- 
fined by 



Up to a constant matrix coming from a martingale term (due to the noise) and thus 
irrelevant for us (see [2], [7]), the coefficient £! is defined by the Green-Kubo formula 



The signature of the superdiffusive behavior of the system is seen in the divergence of 
$), i.e. in a slow decay of the current-current correlation function. To study the latter we 
introduce its Laplace transform 



which is well defined for any z > 0. 

Our third theorem is the following lower bound on &(y,£). Observe that &(y,z) is a 
square matrix of size 2 whose entry (i,j) is denoted by J^ij. 

Theorem 3. Fix y > 0. There exists a positive constant c := e(y) > such that 




(11) 




(12) 




&i,i(Y,z)>cz-V 4 



and 



^ij(y,z) = 0, (i,j) ? (i,i). 

Moreover, there exists a positive constant C := C(y) such that for any z > 0, 

r^i,i(i ) z/r)<^i(y,z)<c^ |1 (i, z //), 



(13) 
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The last part of the theorem is easy to prove but has an important consequence. In [7] 
numerical simulations are performed to detect the anomalous diffusion of energy. Since 
it is difficult to estimate numerically the time autocorrelation functions of the currents be- 
cause of their expected long-time tails, a more tenable approach consists in studying a non 
equilibrium system in its steady state, i.e. considering a finite system in contact with two 
thermostats which fix the value of the energy at the boundaries. Then we estimate the 
dependence of the energy transport coefficient k(N) with the system size N. The latter is 
defined as N times the average energy current. It turns out that k(N) ~ N s with a parame- 
ter 5 := 8(y) > increasing with the noise intensity y (except for the singular value 5 = 1 
when 7=0 which is a manifestation of the ballistic behavior of the Kac-van Moerbecke 
system). This result is very surprising since the more stochasticity in the model is intro- 
duced, the less the system is diffusive. The same has been observed for other anharmonic 
potentials in [7] and also for the Toda lattice perturbed by an energy-momentum conser- 
vative noise ([19]). It has been argued in [19] that this may be explained by the fact that 
some diffusive phenomena due to non-linearities, like localized breathers, are destroyed by 
the noise. In [3] simulations have been performed directly with the Green-Kubo formula 
for other standard anharmonic chains with the same conclusion: current autocorrelation 
function decreases slower when the noise intensity increases. If all these numerical simu- 
lations reproduce correctly the real behavior of the models investigated, they dismiss the 
theories which pretend that some universality holds, e.g. [29]. It is therefore very impor- 
tant to decide if the phenomena numerically observed are correct or not. The inequality 
(13) shows that the time decay of the current autocorrelation function is independent of y 
(up to possible slowly varying functions corrections, i.e. in a Tauberian sense). Therefore 
the numerical simulations do not reflect the correct behavior of the system 3 . Of course 
we could imagine that the exponent 8 appearing in the transport coefficient defined by the 
non-equilibrium stationary state depends on the intensity of the noise and not on the au- 
tocorrelation function but this is very improbable. The last part of our theorem is in fact 
valid for all the models cited above. It applies in particular to the framework studied in 
[3] and shows that the numerical observations of that paper are not consistent with the real 
behavior of the system. 



To study the energy-volume fluctuation field 3( n ' a , we introduce the following change 
of variables S, x = e~ br,x , for each x € Z. Then, the previous Markovian system (tj(*))<>0 
defines a new Markovian system [E, (f ))r>o with state space (0, +°°) z whose generator j£f 
is equal to b 2 s/ + yS?, where for local differentiable functions / : (0, — > R we have 
that 



Observe that the energy and volume conservation laws correspond, for the process 
(§(/))t>0, to the conservation of the two following quantities Y,xeZ%x> L.vezlog^- The 
corresponding microscopic currents are defined by the conservation law equations: 



4. A CHANGE OF VARIABLES 



and 




'It would be very interesting to understand why the numerical simulations are so sensitive to the noise. 
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where 

7^+i ) = -b%Ui - rv^, f XtX+1 (§) = -fe 2 (^ + <?, + 0- rv(iog^). 

We have the following relations between the microscopic currents 

!xs+i(n) =i^+i(^)-^+l(4). /U+iC 1 ?) = ( 14 > 

If 77 is distributed according to jXp ^ then £ defined by <i; x = e ^ is distributed accord- 
ing to the probability measure ^ on (0, +°°) z given by 

Vp^) = nZ" 1 ()3,A)l { ^. > o } e X p{-j3^ + Alog(^)} 

xez 

with Z(j3 , A) the partition function and 

J3=J3, A = -l+A//?. (15) 

Remark that Vp x is nothing but a product probability measure whose marginal follows 
a Gamma distribution Yx+irf- 1 w ^ parameter (A + 1 , y3 1 ) . In particular, we have Z := 
Z(j3,A) = j3~^ +1 'r(A + 1), where r is the usual Gamma function. 

Thus, the process (f ))r>o has a family of translation invariant measures Vp x parame- 
terized by the chemical potentials (p\A) <G (0,+°°) x (— 1,+°°). 

Let P Vj3 x be the probability measure on D(R+, (0, +°°) z ) induced by (| (t)) t >o and Vp x 
and let E V/j x denote the expectation with respect to P Vjg k . 

Let {-)vp X denote the average with respect to Vp x- The averages p := p(j3,A),0 := 
9 (j3 , A ) of the conserved quantities for (£, (?))?>() at equilibrium under ^ are defined by 
P = (^x)vfl A and = (log(^)) Vj3 A . By a direct computation we get 

p = l + e-bv = ^-, = -Z>v = y(A + l)-log(j3) (16) 

where is the usual digamma function, i.e. \j/ = T' /T. 

It is understood, here and in the whole paper, that (p\A) are related to (p\A) through 
(15). We will use the following compact notation, for each x € Z, 



log(fe) 



[ | 



a=( i ,1 i. 



-l/b 

The co variance matrix % : = x(P, A) of (Oq under ^ is given by 

d|log(Z) -^ iA log(Z) \ / (A + l)/3- 2 >3 1 



* V -^iog(z) ^ 2 iog(z) J { -p- 1 (lo g r)"(A + l 

Thus, the covariance matrix % of (Bo under Vp x is related to the co variance matrix % of 
<»o under fip^ , by 

Z=A^A r . (18) 
A simple computation shows that (j x . x + \ ) Vj3 A = —b 2 p 2 and O'^+i)^ = —2b 2 p. 
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The hydrodynamical equations for the process (4 (f ))r>o are given by 

\d t p-b 2 d q {p 2 )=Q 

\d,e-2b 2 d q p = 

and can be written in the compact form d t X + d c/ 3(%) = with 



(19) 



The differential matrix of 3 is given by 

(va)(x) = 



-2b z p 
-2b 2 



As above, let (T t + ) t >o (resp. (7J )t>o) denote the semigroup on 5(K) x S(R) generated by 
the linearized system 

d,e+M T d q e = 0, {resp.d t e-M T d q e = 0). (20) 

where 

M:=M(p,0) = (V3)(£»), 

p and are given by (16). We omit the dependence of these semigroups on (p,9) for 
lightness of the notations. 

We remark that the linearized system of (19) around the constant profiles (p , 0) is given 
by the first equation on the left hand side of (20). It is easy to show that M = AM AT 1 and 
A T f- = T,-A T . 

5. Orthogonal decomposition 

Observe that vp x is a product of Gamma distributions. Let us recall that the Gamma 
distribution Jajc with parameter (a,k) is the probability distribution on (0, +°°) absolutely 
continuous w.r.t. the Lebesgue measure with density f a ^ given by 

f a , k (q) ={k a T{a)y\ a - l e-ll\ q > 0. (21) 

Thus, we have V W (rf|) = Uxez (fi+i,p-i&)d$ c ) = FL e z (PA+l,l(P&)d&). The 

generalized Laguerre polynomials (H„) n >o form an orthogonal basis of the space L 2 (7a+i,i)- 
They satisfy the following equations: 

flW = 1, 

(q^L+ (A + 1 - q )^- + n)Hi X) = 0, 
V aq dq / 

(n + 1 (q) = (2« + 1 + A - q)H™ (q) - (n + A )H^\ (q) 
and the normalization condition 

r(A+» + i) l 



(22) 



jf " (//<%)) fx+X,l{q)dq = 



r(A + i) «!' 
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In particular, we have 



(A) (2 + A)(l+A) q 
r 2 W = 2 ( A + 2 )<? + ^ 



2 (23) 



Let E be the set composed of configurations a = {<J x ) x ei G N z such that a x ^ only 
for a finite number of x. The number £ ve z °* is called the size of a and is denoted by |<r|. 
Let E„ = {ff £ L ; |cr| = «}. On the set of n-tuples x := (x\, . . . ,x n ) of Z", we introduce the 
equivalence relation x ~ y if there exists a permutation p on {1 , . . . , n} such that x p ij\ = y, 
for all i <G {1, . . . ,«}. The class of x for the relation ~ is denoted by [x]. Then the set of 
configurations of E„ can be identified with the class of n-tuples classes for ~. 

To any (7 <G E, we associate the polynomial function given by 

Then, the family e} forms an orthogonal basis of L 2 (v M ) such that 

H'W/dVv = ^IL ^1 |=^), (24) 
where 

on „ T(A + <7 r + 1 ) 1 

:=ILs V(A + 1) ^ (25) 

and 5 denotes the Kronecker function, so that 8 a=a i = 1 if a = o' , otherwise it is equal to 
zero. 

The functions F such that F(a) = as soon as |<y| ^ n are called degree « functions. 
A function F:E^1 such that F[a) = if (7 ^ E„ is also called a degree n function. 
Any function F : E„ —> M is nothing but a symmetric function F : IP — > M through the 
identification of a with [x]. We denote by (-, •) the scalar product on ©L 2 (E„), each E„ 
being equipped with the counting measure. Hence, if F, G : E — > R, we have 

(F,G) = £ £ F n (a)G n (a) = £ I £ F n (x)G„(x), 

«>0ctgI„ n>0 n- xGZ" 

with F n ,G n the restrictions of F, G to E„. Since (/3,A) are fixed through the paper we 
denote /Z^" 1 by // CT . 

If a local function / e L 2 (v | g ^) is written in the form /(<!;) = LctgI^( (J )^(^) men 
we have 

= I (2tF)(cr)H CT (§), = £ (6F)(a)H ff (§) 

with 

(6F)((T)=£(F((T-+ 1 )-F(CT)), 

*ez 

where o' Jr - v + 1 is obtained from a by exchanging the occupation numbers a x and G x +\. 
Let us now compute the operator 21. We have 



.teZ 



II 



By the definition of H c and by the second equality in (22), it follows that 

agZ 

where a — 8 X is the configuration where a particle has been deleted at site x (if there was 
no particle on site x, then a — 8 X = a). 

Now, noticing that the fourth equality in (22) can be written as 

(/?<?) = (2n + 1 + A - (n + k)H^\(Pq) {n+l)H^{f$q) 
and performing some change of variables, we have that 

(J*» ff )(§)= £ + A)(<* + l)*a+$-«<(5) 

x,yeZ 
\x-y\ = \ 

xel 

+ E + 1 ) - <*v- 1 ) ff CT +5« )• 

agz 

Here, a(z) = — 1 if z = — 1 , a(z) = 1 if z = 1 and otherwise. It follows that 

2l = 2l + 2U+2l + 

with 

(21oF)(ct) = - £ fl (y-x)cj v ((j y + l+A)F((j + 5 v -5,), 

A.yGZ 

|a-v| = 1 

(B+F)(a) = -£ cr^cTv+i - p^i)F(a - 

agZ 

(SUF)(a) = £ (<t t - 1 + A)(a v+1 - a^i)F(ff + &). 

xel 

Observe that if F vanishes outside of £„ then %±F vanishes outside of E„ T i and 2lo 
vanishes outside of £„. 

The Dirichlet form &{f) of a local function / <G L> 2 (Vj3 a) is defined by 

^(/) = (/,(-^/))v w =^I / (/(^ ,v,+1 )-/(^)) 2 v /5 ,a(^). 

agZ 7 

Here and below (-,-)v^ A denotes the inner product of h 2 (vpx)- 
If / has the decomposition / = LaGE^X* 7 )^ then 

9{f) = W I ^((T)(F((T- +1 )-F((T)) 2 . (26) 

aGZ cgL 

Let A+ = { (x,y) G Z 2 ; y > x + 1 }, A_ = { (x,y) G Z 2 ; y < x - l} and A = {(*,*) ; x G Z}. 
We denote by Di the Dirichlet form of a symmetric simple one dimensional random walk, 
i.e. 




where F : Z -> K is such that LagZ-F 2 ^) < +oo. 
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We denote by Di the Dirichlet form of a symmetric simple random walk on 1? where 
jumps from A± to Ao and from Ao to A± have been suppressed and jumps from (x,x) G Ao 
to (x± l,x± 1) G Ao have been added, i.e. 

W = ^£ I (F(x + e)-F(x)) 2 + W(F(x±(l,l))-F(x)) 2 , 

|e| = lxeA ± ,x+eeA ± xeA 

where F : Z 2 — > R is a symmetric function such that £ x gZ 2 F 2 (x) < +°°- 

Lemma 1. Let f = YJn=\ £crei„ F n (<5)H a be a local centered function of degree less or 
equal to 2. There exists a positive constant C := C(A), independent of f, such that 

C- l [D l (F l )+B2(F 2 )}<&(f)<C[J} l (F l )+B2(F 2 )]. 

Proof. Observe that 

• If o- G Zi, then = (A + l). 

• If a G E 2 , a = 8 X + Sy, x^y, then W{a) = (A + l) 2 ; if a G E 2 , a = 2<5 r , then 
^(ff) = [(A+2)(A + l)]/2. 

This follows from the relation T(z + 1) = zT(z). Then, by using (26) and the identification 
of functions F : E„ — > R of degree « with their representations as symmetric functions on 
Z", the claim follows. □ 

6. Triviality of the fluctuations 

In this section we prove Theorems 1 and 2 and the three Corollaries above. Since 
Vb{r\ x ) — 1 = %x, — l°g(£vX r \x = — b 1 l°g(£x)i the problem is reduced to study the fluctu- 
ation field of the conserved quantities for the process {t, (f ))f>o at equilibrium under the 
probability measure Vp x. The latter is defined by 

>«( G ) = -L £ G (jt/n) • (co x (tn l+a ) - co) , (27) 
where G is a test function belonging to S(R) x S(R). Recalling (17) we have 

^».«( G ) = -L £ (A 7 G)(*/«) • (flfc(fn 1+B ) - co) = Jf '"(A r G). (28) 

v n xez 

By the relation M = AMA -1 , we are able to translate any result about the convergence 
of 3f."' into a corresponding result for . 

6.1. The hyperbolic scaling. For any local function g := g(£,) we define the projection 
^p,e g °f g on the fields of the conserved quantities by 

{& p ,eg) (§) = *(£)- Kp, 0) - ( V D (P, 0) ■ (coo - co) 
where g(p,9) = (g)v pl and Vf is the gradient of the function g . 
We have that 

Proposition 1 (Boltzmann-Gibbs principle I). For every local function g G L 2 (vpx), far 
every H G 5(R) x 5(R) and every * > 0, 

\ 2" 



lim E v „ . 



^H(x/«).[T r ^ p ,^(^( OT ))] <fe 
o V« 



= 0, 



•xez / 
where for a local function g and | G (0, +°°) z , T v g(<i; ) := g(t x % ) fl "^ Tv^ (y) := 4 (* +>')• 
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Since we prove a refined version of this proposition we omit its proof. As a consequence 
of last result, we get that the fluctuation field (i2:"' )„>i converges in law (in the sense of 
Theorem 1) to the solution of the equation at the right hand side of (20). Theorem 1 is a 
simple consequence of this fact. 

In order to prove Corollary 1 and Corollary 2 we follow the same argument as in the 
proof of Theorem 4. 2 of [ 1 8] . In order to compute the limiting variance for the energy flux, 
we do the following: 

E e [Z r %] = EqIWJ?,}) - %{H l )}{W s {Hl) - %(H*)}] 

= 1brE q [9((G\)91(G\)- & t {G\)%{G\) - W S {G\)%{G\) + %{G\)%(G))] 



where 



Ho(x) \ n i f ,_(G t {x) 



Hq is the Heaviside function, namely Hq(x) :~ l{ x >o} and G((x) := (1 — -*/^)l{o<x<f}- By 
(28), E Q [^ t (H)%(G)] := (Tf (A T H) ■ ^A r G). By the definition of (Tf) t > we have for 
G\,G 2 test functions in ; 



T ( G 1 (x)\ = l -[G 2 {x-2b 2 pt)-G 2 (x))+G l (x-2b 2 pt) \ 
' \ °2( x ) J \ G 2 (x) J ' 

As a consequence we obtain that 

x lim / (G\{x)G\{x) - G\{x)G t {x) - G s e {x)G e (x) + G e {x)G e (x))dx, 



Se[2*zj]= l- 



where for t > 0, G[(x) := Ge(x — 2b 2 pt). Now, using (15) and (16) the proof ends. Anal- 
ogously, by taking an approximation vector function as (0,//o) we get the result for the 
volume flux. 

6.2. The longer time scale. Since in the hyperbolic time scale the initial fluctuations for 
the field 3f."' a are transported by the linearized system on the right hand side of (20), we 
redefine the fluctuation field S h - a , a > 0, on G G S(R) x S(R), by 



J° ' (G) = Sfp \T t „aG 



For t, G (0, +oo) , let J Xr x+i ) be the microscopic current vector 
By Dynkin's formula, 

.<'"(G) = if- a (G) - J*"(G) - J' \n x+a J? (if' a (G)) +d s ^'- a (G)} ds 

is a martingale with quadratic variation given by 

(JZ n ' a ) t = f Q nl+a ^ (^'"(G)) 2 -2« 1+a (^''"(G)) Jgf (if' a (G)) ds. 

A simple computation shows that E V)3 x [(^#" ,a )r] vanishes as n goes to +°° for a < 1. This 
is equivalent to saying that the martingale „#"' a vanishes as n goes to +°° in L 2 (P V ^ A ), 
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for a < 1. Observe that, by definition of (r, + ) r >o, we have 



V n xeZ 



sn l+a )-a>) 



£ [d q {T+ a G){x/n)\ .M(0) x (,« 1+a )-C0). 



On the other hand, the first term in the integral part of the martingale ./#"' a (G) is equal to 

- &) ( - ) ) • (■/,,,+! (<^« 1+a )) - Ev M [/^i] 



Performing a Taylor expansion, we can replace this term, up to a term vanishing as n goes 
to +00 in L 2 (P V ),by 



>/3.A 

,a 



((d q T s + n aG)(x/n)) ■ (/,., + i - (Vn^A 



Thus, in order to show that 



3f?' a (G)-2%' a (G) 



(29) 



it remains to show that 



lim E v 



li.A 



in Jo 



dsZ(d q T+ a G)(x/n)-® x (!;( S n 1+a )) 



= 



where for % e (0, +°°) z 

= (§) - (/^i)v w -M((0^- ft)). 

Observe that in this formula, M := M(p, 9) is the differential w.r.t. (p, 0) of the function 
(J x ,. v +i)vp A as computed below (20). A simple computation shows that for t, e (0, +°°) z 



-* 2 (&H-i-p)(&-p)-(y+ft 2 p)vfe 

-V^ + ylogO 



The discrete gradient terms appearing in the previous expression, permit to perform 
another discrete integration by parts and the resulting terms vanish in L 2 (P V/j x ) as n goes 
to +00, for a < 1 . Using the smoothness of the function G, we see that it only remains to 
show the following 



Theorem 4 (Boltzmann-Gibbs principle II). Fix a < 1 /3 and let (p : M + x 

that <p(-,x) € S(M). For every t > 



: be such 



lim E, 



/J .A 



t „a 



+ 1 ( OT 1+a )- P )^ 



Proof. In the following, C,Co,Ci , . . . denote constants independent of n whose values can 
change from line to line. Let f s {^) be the function Y,xcz ( P( s : x / n )Hs x +s x+l (§) ■ We have 
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the following upper bound 




In the first inequality above we used the result of [25] applied to this setting. We notice 
that since our test functions depend on time, the result of [25] has to be modified as written 
here. To prove the last result one can simply adapt the proof of Lemma 4.3 of [10] to this 
case. We also notice that for our purposes, we need to prove the result for <p equal to the 
first component of the column vector d q T^ a G. 

In order to simplify notations, let us define e = 1 /sn +a . 

We denote by 1$, the set of configurations a of £2 such that a = 28 x , x <S Z, and the 
complementary set of m ^2> i- e - the set of configurations a € £2 such that a = 8 X + S y , 
y^ieZ. Observe that f sn a is a function of degree 2 with a decomposition in the form 
fsn a = Laei.2^.sn a {^)Ha which satisfies <S> sn a(o) = if a <E £3. We have that (see e.g. 
[25]) 

(fsn°>,[ — 1+^-7^] /««) =SUp{2(/ OT «,g) v s k -e(g,g) v - 
x \ *" / ' V /3,A g J 

where the supremum is taken over local functions g £ L 2 (v^a)- Decompose g appearing 
in this variational formula as g = £ CT G(o)H a . Recall that {H c ; a G £} are orthogonal, 
that the function f sn a is a degree 2 function such that <& sn a (a) = for any a ^ £j and 
formula (26) for the Dirichlet form ^(g). Thus, we can restrict this supremum over degree 
2 functions g such that G(<j) = if a £ £?. Then, by Lemma 1, we have 




where C := C(A,y), A± = {{x,y) £ Z 2 ;i ^ y} and as usual we identify the functions 
defined on £„ with symmetric functions defined on II 1 . 

To any symmetric function G defined on the set A±, we denote by G its extension to I? 
defined by 

G(x,y) = G(x,y) if x ^ y, G(x,x) = -£ G{(x,x)+e). 

\e\=l 
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It is trivial that 

£ & 2 (x,y)<C £ G 2 (x,y), 

(x,y)eZ 2 (x,y)eZ 2 



£ £ (G(( X ,y)+e)-G( X ,y)) 2 <CZ £ (G((x,y) + e) - G(x,y) 



\e\=l (x,y)eZ 2 

Thus, we have 



|e| = l to)6A± 
(.tjO+eGA^ 



)1A 



<C sup^ £ 4>. v „a(x,y)G(x,y)-Ci£ £ G 2 (jc,y) 



(,v,y)eZ 2 



-C2E £ (G((x,.y)+e)-G(^) 

|e|=l (x,>)6Z 2 

where the supremum is now taken over all symmetric local functions G : 1? —> K. By using 
Fourier transforms, it is easy to show that this supremum is equal to 



dk 



Cor |4>,„, a (k)| 2 
4 i[o,i] 2 Ci£ + 4C 2 I 2 = i sin 2 (7r^ ; ) 

where the Fourier transform <t> s „a of <t> v „a is given by 

4> s „«(k)= £ ® s „ a (x,y)e 2i "^ x+k2 y\ k=(fe 1 ,fe 2 )G[0,l] 2 . 



By definition of f sn a, we have 4> ln a(jc,y) = -\q>{sn a ,x/n) + (p(sn a ,y/n) ) if |jc — y\ = 1 
and otherwise. Consequently, we have 



-1 



fsn a ! 



,1+a 



fsn" ) < S 



E* 6 z (p{sn a ,x/n)e 2illx( - kl+k ^ 



Co 


( < 


~ 16j 


[0,1] 


" 16j 


( . 

[0,1] 


C 

s 7? 


/ 

i[0,l] 



v,j.a 167[ai] 2 Ci£ + 4C 2 £ 2 =1 sin 2 (7rfc,) 

2 



-dk 



£ p(an°,Jc/»)« 2<,DSP 
reZ 

£ 9(i« a ,x/n)e 2top 

a-gZ 

Observe now that 



[0,1] Ci£+4C 2 sin 2 (7rfc 1 ) +4C 2 sin 2 (n;(/7-fc 1 )) 
dk\ 

[0,1] C i £+4C 2 sin 2 (7rA: 1 ) 



dp. 



'[0,1] 



(-t,y)GZ 2 



= £ <p(x/n)<p(y/n) / ^^"dp = £ <p 2 (x/«) < Cn. 



[0,1] 
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Putting everything together, we get that 



E, 



Ctn 2a - 1 f n 
~ ~~tA+a~ Jo Ve ' 

Since e := 1 /sn l+cc last expression vanishes as n goes to +°o, if a < 1/3. □ 

Now, in order to prove Corollary 3 we follow the same arguments as in the proof of 
Proposition 9.3 of [18] and we proceed as follows. We use the change of variables to 
defined the energy (resp. volume) flux through the time dependent bond uf a («) during the 
time interval [0 7 tn l+a ] as: 



<Ho (f):= E {^(m 1+a )-^(0)} 



y>u,' (n) 



resp. r% (n) (r):= £ {log&(m 1+a )) - log(§,(0))} 



y>u t ' (n) 



By (14) we have that 



g* a At) := $\a. At) - f'l.a. At), f\a. At) 

Then, from (29) applied to the vector function G\ we obtain that 



-\fla. At) 

b u i W 



lim E v 



^{i\ a , At)-E VR ,[£", ai At)}} 



= 0. 



On the other hand, applying (29) to the vector function 

±G e (x-uf a (n)) 
—Ge(x— uf a {n)) 



we obtain that 



lim E v , 



/J.A 



Now, Corollary 3 follows easily from the previous results. 

7. DlFFUSIVITY 

In this section we prove Theorem 3. Our proof is based on the resolvent methods intro- 
duced in [4, 22] and developed in few other contexts (e.g. [5, 24, 28]). Some differences 
with these previous works is the presence of two and not only one conserved quantity and 
the degeneracy of the symmetric part of the generator. 

We fix p > 0, 6 G K and denote by p\A the chemical potentials given by (16). Let also 
(p\ A) be given in terms of (jS , A) by (15). 

Recall the definition of J XyX +\ given in (1 1). We introduce the normalized currents j AVC +i, 
i'xx+\ ar, d J-v+i corresponding to the process {£, (t)) t >o, which are defined by 

j.r r r+l(4) = J^+l(|)- {jx^c+l)v PtX -^p(;^+l)v M (^.v-P) -<5e(A,.v+l)v j3 A (log(§c) - 0), 

= y^+i(4)- (iw+i)v f ,, -^p(A r x+i)v fiA {£,x-p)-deU'x,x+\)vp A (iog(<y - e), 

J*,*+ 1 ) = 1 (4 ) - 1 (I ) ■ 



= 0. 
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Since (j x , x+ i) Vpx = -b 2 p 2 and {f x ^ +l ) Vfj x = -2b 2 p, we get 

j„ + i(|) = -b 2 (^ x -p)(Ui -P) - (y+* 2 p)V(£,) 

(30) 

j;* + i(^) = -v(^, + yiog(£ T )). 

For any local compactly supported functions f,g : (0,+°°) z — > K we define the inner 
product < /,g >:=< f,g > |Sa of / and g by 

((^)v M -(/)v^(g)v M ) 




where the third equality follows from the invariance of Vp x by the shift. Observe also that 
the first sum on Z is in fact a finite sum since / and g are assumed to be local functions. 
We denote by ,3% the completed space. Observe that any constant and gradient functions 
are equal to in ,3%q. 

By (30), the normalized current associated to the volume is a gradient and this shows 
that ^,j(y,z) = if ^ (1,1). By the definition of J x . x+ i and by (14), we are only 
interested in the behavior, as z — > 0, of 

£(z) =<Jo,i,(z--2T 1 Jo,i »= / e ~~ J <Jo,i(O.Jo,i(0)» dt. 

Jo 

Since gradient functions are equal to in J^, this is equivalent to estimate 
£(z) = b 4 « Wb.i , (z - %)~ X W ,i » 

where W VjV is the local function W V;V = — p)(| v — p)- 

In this section we prove there exists a constant C > such that 

< W ,i,(z-JS?)- 1 Wb,i >> Cz~ 1/4 . (31) 

But before proving (31) let us show (13) which is a direct consequence of the following 
lemma. 

Lemma 2. For any y > 0, f/iere ex;'ifi a constant C := C(y) smc/i f/zaf 

<w ,i, (z/y-fcV-^ , )- 1 w oa >< C<W ,i,(z-/7V-7^)- 1 W oa > 

one/ 

<Wo.i,(z-^V-7J^)- 1 W ,i >< C<Wb,i, {z/r-b 2 £/-,Y)- 1 W( ) A > . 

Proof. Assume y > 1 the case y < 1 being similar. By Lemma 2. 1 of [4] we have the 
variational formula for <C Woj, (z — Jtf)~ l Wo.i where _Sf = b 2 si ' — yS?, given by 

sup {2 « Wo.! , / » - « /, (z - yy)f » -& 4 « rff, (z - y^)- W/ »} , 
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where the supremum is carried over functions / belonging to the domain of the generator 
S£ or equivalently to a dense subspace included in this domain, say the space of smooth 
local compactly supported functions. We have that 

sup {2 « W Q . u f » - « /, (z - yy)f » -b 4 « rff, (z - yJ?T V/ »} 
/ 

= sup{2«W , 1 ,/» -y</, (z/y-^)/> -&V 1 «<f,(z/y-^T 1 ^/»} 
> sup { 2 « Wo. i , / » - r « /, (z / y - ^)/ » r « *«7, (z/y - " 1 <f » } 
=su P {27- 1 / 2 « w,).! ,/ » - « /, (z/y- ^)/ » -fe 4 « (z/y- V/ »} 

where the penultimate inequality comes from y > 1 and the ultimate is obtained by the 
change of / into y~ l / 2 f- The last term is equal to 

y" 1 «W ,i, (z/y- ftV- ^) _1 W 0i i > 

and this proves the first inequality of the lemma. 
For the second one we proceed similarly: 

sup {2 « W Q . u f » - « /, (z - y^)/ » -fo 4 « (z - y^)"W/ »} 
/ 

=su P {2 « w ,i ,/ » -r« /, (z/y- -y)f » -^V 1 « (z/y - -*TW »} 
/ 

< sup { 2 « Wo. i , / » - r 1 < /, (z/ 7 - ^)/ » -ft V 1 « stf, (z/ 7 - y) ~ ' *7 » } 

/ 

[ly 1 ' 2 « Wo,i,/ >-«:/, {z/y-y)f » -ft 4 « (z/y- ^)"V/ »} 



= sup 

/ 

=7<Wb,i, (z/y-fcV-^) _1 W ,i > 



□ 



Recall the orthogonal decomposition obtained in Section 5. Let / = £^^((7)//^ and 
g = L <7 G(o')// CT be two centered local functions (i.e. F(0) = G(0) = 0, where denotes 
the empty configuration of £). The shifted configuration a by z £ Z is denoted by T : a. We 
identify F ni G„, the restrictions of F, G to £,,, to symmetric functions on Z". By (24) we 
have that 

«/,g»= E E F(T z C7)G(0)>r(C7) J 

where #^ was defined in (25). 

With some abuse of notations, we denote by <C F, G ^> the scalar product defined by 

<F,G>= E E F (T z C7)G(C7)>r(C7). 
zeZ ctgE 

We also introduce the inner product <C •, • ^>f re e defined by 
<F,G> ftee = E E *"(ijff)G(a). 

Since the function W is invariant by the shift, we have a very simple relation between 
these two inner products: 

<F,G> = < yr 1/2 F,#- 1/2 G> fr ee . (32) 
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On the set E« we introduce the equivalence relation * defined by (7 * o' if and only if 
there exists «gZ such that t u a = o' . Let E* = £„/* be the set of classes for this relation 
and E* = U„>iE*. We can rewrite the scalar product ■, ■ 3>free as 

<F,G» fre e= £ F(o)G(a). 

aeZ* 

Here F is defined by F(ct) = (Y, ye z TyF) (a) where a is any element of a. The function 
W being invariant by the shift, it makes sense to define W(o) by W(o), ff £ o, a € E*. 
Then, we have 

«F,G>= £ # / (a)F(a)G(a). 

Lemma 3. There exist constants C :— C(n,X),c := c(«, A) such that for any local function 
F : E„ — > K of degree n and any positive real z > 

CT 1 < F,F >f ree < < F,F > < C < F,F > free • 

(2) 

CT 1 < F, -&F > free <«F, -6F > < C < F, -6F >f ree ■ 

(3) 

< F, ( z - yer'F >=< #- 1/2 F , (z - 76)- 1 ^ 1/2 F > free . 

Proof. Recall the definition of W from (25). Thus, W is bounded from above (resp. from 
bellow) by a constant C(n,X) (resp. C~ l («, A)) independent of a G E„. This is enough to 
conclude (I). In order to prove (2), it is enough to use (32) and the fact that for any local 
function F:E->lwe have that &(W l/2 F) = W l/2 &F. Finally, for a local function F of 
degree n, we have by (32) and the fact that 

<F,(z-7©)- 1 F>= sup {2<F,G>-<G,(z-7©)G>}, 

G of degree n 

the following equality 

< f, (z-y©)-^ >=< 5r'/ 2 F (z- re)- 1 ^ 1 /^ > free , 

which proves (3). 

□ 

Our goal is to get a lower bound for <C Wo,i, (z — -Sf ) _1 Wb,i 3> which by Lemma 2.1 of 
[4] can be rewritten in the variational form 

sup {2 « w ,i,f >-«:/, (z- r^)/ » < *7, (z - y^r 1 ^/ »} • 

/ 

Any element a G E* can be identified with an element of 1 through the application 
which associates to ((X\,...,(X n -\) G N the class of the configuration a = 5q + 8 a] + 
■ ■ ■ + 8ai+...+a„_ l - 

Observe also that 6 is a self-adjoint operator with respect to <C -, ■ ^> and with respect 

to < •,• >free- 

We restrict the previous supremum over degree 2 functions / = Y,t X y)eZ 2 ^i\ x ^y\)^[x,y] ■ 
In order to keep notation simple, whenever we identify a configuration a G E„ with [x] G Z" 
we will simply write F(x), instead of F([x]). 
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Up to some irrelevant multiplicative constant, a lower bound is given by 

sup {2F(0,l)-||F||^-^||2(_F||^, c -^ 4 p + F|| 2 _ 1 , ; -^ 4 poF|| 2 _ 1 , z } 

Fof degree 2 

where ||F|| 2 j , =<§; F, (z — y&) ±l F 3>. We also introduce the corresponding H±i , -norms 
associated to < • > free : ||F||| Ufree =< F, (z - y&) ±l F > fre e, for F : E -> R. 

By Lemma 3, there exists a constant C such that this lower bound is bounded from 
bellow by 

sup {2F(0,1) -C||F|| 2 +lAfree 

Fof degree 2 

Let us first show that if F is of degree 2 then the contributions given by ||#' 1 ' 2 2l_F|| 2 _ 1 f 

and || W 1 ^ 2 ^IqF\\ 2 _ j Tfree are equal to zero. 

The function W is constant equal to (A + 1) on Ei so that W X l 2 %-F = \/A + 12t_F. It 
is easy to check that the degree one function 2l_F satisfies 

(2l_F)(w) = (A-1)(f(w-1,k)-F(k,k+1)). 

For any degree one function G, we have 

«2l_F,G> free = £ G(M + y)(A-l)(F(u-l,u)-F(M,ii + l))=0 

by a telescopic sum argument. This shows that 2l_F is equal to zero in the Hilbert space 
generated by <C •, • >f ree - 

Recall that if F is a degree 2 function, i.e. a symmetric function on Z 2 , then F is 
identified with a function F defined on N by 

F{a) = £ F(u,u + a) 

uez 

and as a consequence, for F and G degree 2 functions it holds that 

<F,G> ftee = £F(a)G(a). (33) 

aeN 

Observe that (21oF)(k,v) is equal to 

'2(1 + A)^F(m-1,m)-F(m,h + 1)Y ifw = v, 

(1+A)^F(m-1,m + 1)-F(m,m + 2)) + (2 + A)(f(m,i<)-F(m + 1,m+1)), 

if (m, v) = (u,u + 1), 
(1+A)(f(m-1,v)-F(k+1,v)+F(k,v-1)-F(m,v + 1)Y if |k-v|>2 

and 

W(u,u)= ( A + 1 K A + 2 ) ; #-( M;V ) = (A + l) 2 fort/^v. (34) 
It is then easy to show that 

W l / 2 {^ F)(a) = 
for any a € N. Putting together the previous result and (33) it follows that: 
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= £ ^ l/2 (%oF)(a) (A-7e)- 1 ['r 1 /2(2l 0j F)]( a ) =0 . 

ceeN 

Lemma 4. There exists a positive constant C such that for every symmetric function F of 
degree 2, ifF(a) = E z eZ^(z,Z+ 05 ), f/ie« 

C" 1 £ (P(y)-F(x)) 2 <<F,-6F»fiee < C £ (f(y)-F(jc) 

|.t— .y| = l ;y|=l 

Proof. This follows easily from the following equalities together with (33): 
6F(0) = £(SF) (?,>>) 



= £ (F(y + hy + l)-F(y,y)) + (F(y-l,y-l)-F(y,y)) = 0, 

6F(l)=£(<3F)(y,y + l) 

yeZ 

£ (f (y- l,y + 1) -F(y,y + 1)) + £ (^(y,y + 2) -F(y,y + 1)) 



= 2(f(2)-F(1)), 

eF(a) = Y,( &F )(y>y+ a ) 



-- £ ^(y-l,y + a )-F(y,y + a)) + £(F(y + l > y+a)-F(y > y+a) 

yeZ yeZ v 

£ (/?(y,y + a+l)-F(y,y + «)) + £ (F(y,y + a-l)-F(y,y + a 

yeZ yeZ 

= 2fF(a + l)-f , (a)) +2(F(a- 1) -F(a)Y a > 2. 



□ 

To any degree three function G, i.e. a symmetric function G on Z 3 , the function G is 
identified with a function on N 2 : 



G(u,v) = £ G(y,M+y,M + v + _y). 



Since G is symmetric on I? , then G is symmetric on Z 2 . As above, for F and G degree 
3 functions it holds that 

«F,G» fre e= £ F(a,P)G(a,P). (35) 

(aj3)eN 2 

Let D3, acting on the local functions on N 2 , be defined by 

D 3 (G) = £ (G(u + 1 , 0) - G(«, 0)) 2 + £ (G(0, v + 1 ) - 0(0, v)) 2 + (G( 1 , 0) - G(0, 1 )) 2 

W> 1 V> 1 

£ (g(u + 1 , v) - G(u, v)) 2 + (G(u, v + 1) - G(«, v] 



»,r>l 



(36) 
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This is the Dirichlet form of a symmetric nearest neighbors random walk on N 2 where all 
the jumps between {0} x N and N* x W 4 , all the jumps from N x {0} and N* x N* and 
all the jumps from have been suppressed, and a jump between (0, 1) and (1,0) has been 
added. 

Lemma 5. There exists a constant C > such that for any symmetric function G on 1? 

C _1 © 3 (G) <<G,-6G> free < CB 3 (G). 
Proof. We have the following equalities 
6G(0,0) =0, 

6G(0, 1) = 2(<3(0,2) - 0(0, 1)) + (6(1,0) - G(0, 1)) , 

6G(0 ;i S) = 2(g(0,J3 + 1) - G(0,j3)) +2(G(0,j3 - 1) - G(0,j3)) , J3 > 2, 

©G(1,0) = 2(G(2,0)-G(1,0)) + (G(0,1)-G(1,0)) 

6G(a,0)=2(c(a+l)-G(ff,0)) + 2(G(a - 1) - G(a,0)) , a>2, 

6G(a,J3) = (G(a + l,^)-G(a,/3)) + (G(a,/3 + l)-G(a,/3)) 

+ l {a > 2} (c(a - l,p + 1)-G(a,j3)) +l {a > 2} (G(a - I,j3)-G(a,j3)) 

+ l tf > 2} (G(a + l,j3-l)-G(a^))+%> 2} (G(a,j3-l)-G(a^)), a,^ > 1. 

We recognize in these expressions the generator of a symmetric nearest neighbors ran- 
dom walk on N 2 where 

• all the jumps between {0} x N and N* x N*, all the jumps between N x {0} and 
N* x N*, and all the jumps from have been suppressed; 

• a jump between (0, 1) and (1,0) with rate 1 has been added; 

• jumps between (a,j3) and (a± l,/3 =p 1) f° r (f^iP) £ N* x f with rate 1 have 
been added. 

• the non vanishing jumps on N x {0} and on {0} x N have been multiplied by 2. 
This together with (35), imply the lemma. □ 

We choose as a test function F the degree 2 symmetric function F such that 

F(a)=z- l/4 e-^ a - l \ o>l, 

_ ; ~ ' (37) 

F(0)=F(1). 

This function exists since given a function G defined on N we can find a symmetric func- 
tion F defined in 1? such that F = G. For that purpose, take F(x,y) = G(\y—x\)[<j>(x) + 
<p (y)) where the function <p is defined on Z and is such that Y*xez ( x ) = 1/2- Then for any 
a g N, F(a) = Z ueZ F(u,u + a) = G(a)E„ eZ [0(M) + 0(m + a)] = G(a). 

Observe that with this choice, by Lemma 4, 

«F,-©F» free ^z 1 / 4 , ^(l)^^- 1 /^ z ^F 2 (a)^z-'/ 4 . (38) 

It remains to estimate the last contribution given by || W l l 2 G\\ 2 _ l . free where G = 21+ F 
is a degree 3 function. 

4 Here and in the sequel N* := N\{0} 
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Lemma 6. Let G = 2t+F where F is defined by (37). There exists a constant C > such 
that 

ll^ 1/2 G||-l,,free>Cz- 1/4 . 

Proof. For any u, v, w G Z, we have 

G(m,m + 1,m + 2) = / ;, (m,m + 1) — F(h + 1,m + 2), 
G(u,u + l,v) = F{u,v) — F(u + l,v), v > m + 1, 
G(v, m,m + 1) = F(v,u) — F(v,u + 1), v < M, 

G(u,u,u + 1) = 21 F(u,u) — F(u,u + 1) 



G(w, w, u — 1 ) — 2 \F{u — l,u) — F{u,u) ^ 
G(u,v, w)=0 otherwise. 
Let us now compute G(u,v), u,v G N. We get 

G(0,1) = -G(1,0) = 2F(0) - 2F(1), 
G(l,v)=F(v + l)-F(v), v>2, 
G(m,1) =F(u)- F(u+l), m>2, 
G(w,v)=0 otherwise. 



(39) 



By (37) we have that G(0, 1) = G(1,0) = 0. Also notice that G(k,k) = and by (34) we 
have that W 1/2 {u,v) = (1 + A) for u ± v. 

It follows, by Lemma 5, that \\W l l 2 G\\ 2 _ j free is upper bounded by the variational for- 
mula: 



ll^ 1/2 G|| 2 ,, z , fr ee = SUph £ fl(M^ 1/2 (MG( M ,v)-CbB 3 GR) 



sup I 2(1+ A) £ /?(m,v)G(m,v)-C D 3 ( j R) 



R 



(«,v)eN 2 

where the supremum is taken over local functions on N 2 . By (39), we have that 
£ R(u,v)G(u,v) = Y i R(hv)(F(v+l)-F(v))-Y i R(u,l)(F(u+l)-F(u) 

(«,v)GN 2 v>2 u>2 

= £f( V )(r(1, V -1)-R(1,v))-Y,F(u)(r(u-1,1)-R(u,1) 

v>3 w>3 

+ F(2)(r(2,1)-R(1,2) 

= £i?(v)(R(i,v-i)-*(i,v)) -£/?(«)(*(«- l.i) -*(«,i) 

v>2 «>2 

(40) 

We use now the following parametrization of For > 1, v G Z, let us define 

fl(A:,v) =0(&-l,v-&), v>fc, R(u 1 k)^^)(k-\,-u + k),u>k 1 
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where {<j>(k, •) ; k > 0} are functions from Z — > R. We have the following lower bound for 

B 3 (fl): 

3 {R)> £ (r(u + 1,v)-R(u,v)) + (r(u,v+1)-R(u,v)) 

u.v> 1 

which is nothing but the Dirichlet form of a random walk where only jumps connecting 
sites of N* x N* have been conserved. With the choice of the parametrization for R and 
this lower bound, it is not difficult to show there exists a constant C > such that 

n>3(*) >c£ L (<HM+i)-0(M) 2 +(^+M)-<HM) 2 - 

k>0veZ 

The right hand side of the previous inequality is the Dirichlet form of a symmetric simple 
random walk on N x Z. 
By (40), we get 

£ R(u,v)G(u,v) = Y,'l>(0,u)(F(u-l)-F(u) 

(«,v)gn 2 «ez 

where F : Z — > R is defined by = — F(w+2)l{ H >o} — F(l — w)l{ M <_i}. We extend the 
function defined on N x Z to Z 2 by defining <j)(—k,u) = (j>(k,u), k> l,u e Z. Observe 
then that 

nj 3 (/f) > C£ £ (*(*,v+ 1) - </)(A:,v)) 2 + (*(*+l,v) -*(*,v) 

*>0veZ 

= jL I (*(*,v + l)-«(*,v)) + (f(*+l,v)-*(*,v) 

Z <rGZt<eZ 

Consequently we have, for suitable positive constants Ci,C 2 : 



ll^ 1/2 G|| 2 - 1 , are e<ClSUp{2£0(O,«)(F( M -l)-F(«))-C 2 £ (<HU)-<HV)) 2 }. 
tiGZ (u.v)GZ 2 

|u-v|=l 

(41 2 

A standard Fourier computation shows this supremum is of order z~'/ 4 . Indeed, let u 
be the Fourier transform of the function u : IP — > R, defined by 

u(k)= £ e 2i ™\(x), k=(k h ...,k n ), 

xGZ" 

and denote by w*(k) the complex conjugate of w(k). Using the expression of the sum of a 
convergent geometric series, we obtain the following expression for the Fourier transform 
*P(k u k 2 ) of the function (x,y) eZ 2 ^ 8o(y)F(x): 

... ( 1 „-2inh 

which satisfies 

|*P(*i,* 2 )|< . 

Z 3 / 2 + C4Sin"(7Tfci) 

for some positive constants C3,C4. The supremum appearing in (41) is then given by 



1 



C 3 ^ 



l^l,^2)| 



2 



C 9 / = ■ ~ dk\dk2- 

J[oa] 2 z + 4sin"(?z;fci) + 4shr(7z;£ 2 ) 

Then the result follows by a standard study of this integral. □ 
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To obtain (31), by (38) and Lemma 6, it suffices to take a test function in the form aF 
with F given by (37) and a sufficiently small. 

8. Stochastic perturbations of Hamiltonian systems 

In this section we discuss some other possible stochastic perturbations and make some 
connections with the recent models considered in [8]. Let us start with the Hamiltonian 
system (2) with potential V and generator A given by 

A = £(v(iiH-i)-V(n^i)W 

xeZ 

The energy Y*xez ^(l*) and tne volume Y*x&L r {x are conserved by these dynamics. Remark 
that in fact Y,xeZ r l2x and Y,xez "Hix+x are also conserved and that we cannot exclude the case 
that still many others exist. This is the case for example for the exponential interaction for 
which an infinite number of conserved quantities can be explicitly identified. Anyway, we 
are only interested in these two first quantities. The product probability measures x 
defined by 

= Y\Z(j5,X)- l exp{-l5V(ri x ) - A 77.,} drj x , 

xez 

where 

Z(fi,k) = J exp(-PV(r)-Xr)dr. 

are invariant for the infinite dynamics. 

In [7] we proposed to perturb this deterministic dynamics by the Poissonian noise con- 
sidered in this paper and conserving both the energy and the volume. One could also 
consider the " Brownian" noise whose generator S is given by S = Y,xez where 

Y x = (V(ifcH) - V'{r} x -i))d nx + (V'(nx-i)-V'(rix))dr, x+1 + (V'(ifc) - V'fe+i))^., , 
is the vector field tangent to the curve 

{(i?,-i,i?,,i7.x+i)eK 3 ; £ % = 0, £ V(riy) = l}- 

y=x— 1 y=x— 1 

It is easy to see that the process with generator L = A + S conserves the energy and the 
volume and has fip x as invariant measures. A priori, it should be possible to extend our 
result to this system for V of exponential type but the noise S seems to have a quite com- 
plicated expression in the orthogonal basis we used in this paper. The advantage of the 
Poissonian noise is its very simple form. Notice also that the Poissonian noise is a less 
strong perturbation of the Hamiltonian dynamics than the Brownian noise. In this sense 
we are closer to the Hamiltonain dynamics with the Poissonian noise. 

We could also decide to conserve energy and not the volume by adding a suitable pertur- 
bation. The invariant states are then given by , j3 > 0. If V is even, a simple Poissonian 
noise consists to change the sign of r\ x independently on each site x at random exponential 
times. In this case one can prove, as in [6], that the energy diffuses in the sense that the 
Green-Kubo formula converges to a well defined finite value. For a generic V a Brownian 
noise with generator S given by S = Lvgz^v w i m Kx = V' (j\x+i)diu ~ ^'(^-O^v+i ma kes 
the job. 

Consider now the case where we want to add a stochastic perturbation conserving only 
the volume. It does not seem to be easy to define a simple Poissonian noise with such a 
property. A Brownian noise is obtained by the following scheme. Fix j3 > 0, consider the 
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vector field X x = d rjv+1 — <5 7)t which is tangent to the hyperplane {{f] x ^f] x+ \) S ; r\ x + 
rj x+ i = 1} and define the Langevin operator Sp by 

=^l^ 2 +f l(n^ + i)-vu)V, 

where J^ jl = j3 L*ez ^ (^I*) + ^ L*ez It- Observe that 5p depends on j3 but is independent 
of A. The operator 5|3 is a nonpositive self-adjoint operator in h 2 (jip^) for any A and 
Sp (Y, x eZ Ix) — 0- Then, the perturbed volume-conserving model has a generator given 
by 

L v p =A + ySp (42) 

where y > is a parameter fixing the strength of the noise. By construction, the Markov 
process generated by has fip x as a set of invariant probability measures. In fact, using 
the same methods as in [7, 17] one can prove that the only space-time invariant probability 
measures with finite local entropy density are mixtures of the (/X/3.a)a- We can also rewrite 

Ljg as 

L l = I { (l " f ) + # V'ft*) - (l + f ) T'Ofc-l)} ^ 

The microscopic flux y' A ,A+i associated to the volume conservation law is defined by 

Lpiru) = -va-m, h-u = (i + y) Vfa^i) - ( ! - t) y,(7?T) - 

The semi-discrete directed polymer model considered in [8] is, up to an irrelevant scaling 
factor 2, recovered by taking V(rj) = e^, j3 = 1 and y = 2 (see (3.7) in [27]). In [8] the 
authors show that for a particular non stationary initial condition ("wedge"), by developing 
a very nice theory of Macdonald processes, the system belongs to the Kardar-Parisi-Zhang 
universality class ([27]). Unfortunately one can not use their results or their methods to 
derive a more precise picture for the model with exponential interactions considered in this 
paper. For other potentials V the theory developed by Borodin and Corwin in [8] can not 
be adapted but it would be very interesting to see if one can relate the models generated 
by to the semi-discrete directed polymer and deduce some qualitative information from 
the latter. The use of the variational formulas considered in this paper could be the way. 

Appendix A. Existence of the infinite dynamics 

In this section we prove existence of the infinite volume dynamics (| (f ))r>o- We focus 
here on the process t, but the same proof can be carried for the process r\ (or just define 
r\ in terms of ^ by r\ x {t) = —b~ l log^ x (t), ieZ. To simplify notations we will assume 
b=\. 

Since the interaction coming from the deterministic part is quadratic, proving the exis- 
tence of the infinite dynamics is a non trivial task. Nevertheless nice sophisticated tech- 
niques have been introduced by Dobrushin and Fritz in [12]. Here, we follow closely the 
approach of [13] (see also [14, 15]) adapted to our case. By itself, the strategy of the proof 
of existence of solutions is standard: we consider finite subsystems and prove compactness 
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of this family by means of an a priori bound for a quantity E which plays the role of an 
energy density. The obtention of this a priori bound is however non trivial and is the main 
step to get the existence of the dynamics. The aim of this appendix is to show how to 
get such an a priori bound. The a priori bound we derive here for the infinite dynamics is 
also valid for finite subsystems corresponding to a finite set A C 1 with a bound which is 
independent of the size of A. This proves then that the finite subsystems form a compact 
family from which one can extract a subsequence converging to the infinite dynamics. 

We have first to specify the space of allowed configurations £2 C (0, +°°) z . For x £ Z, 
let g(x) = 1 +log(l + |*|) and denote by E(£,H,o), % £ (0,+°°) z , ju G Z, (7 > 0, the 
quantities 

E{S,n,o)= £ (1+2^-log^), 

£(<!;) = sup sup o~ l E(E,,n,o). 

a>g{n) 

The quantity E is called the logarithmic fluctuation of energy and the set £2 is defined 

as 

£2 :={£e(0,+-) z 

The configuration space £2 is equipped with the product topology and with the associated 
Borel structure. It is easy to see that Vp x(Q.) = 1 for any /3 > and A > — 1. 

Let N(f ) = {N XrX+ i (t ) ; x £ Z} be a collection of independent Poison processes of inten- 
sity 7 > 0. The equations of motion corresponding to the generator Jz? read as 

= U^x+i - &-i)dt + V((k - ^)dN x . l>x {t)) , x G Z. (43) 

Let D(R+,R) denote the space of cadlag functions of R+ into R with the Skorohod 
topology and let D = [D(R + ,R)] Z equipped with the product topology and the associated 
Borel field SB. The smallest c-algebra on which all projection restricted to the time interval 
[0, t) are measurable will be denoted by S$ t . Finally, suppose that we are given a probability 
measure P on SB such that our Poisson processes A^+i are realized as components of the 
random element of D. 

Definition 1. A S§ r adapted mapping %(t) := <jj(f,N) mto itself is called a tempered 
solution of (43) with initial configuration ^° G £2 if E,(0) = £ , almost each trajectory 
4(-,N) satisfies the integral form of (43), and the logarithmic energy fluctuation £(£(*)) 
is bounded on finite intervals of time with probability one. 

Theorem 5. For any ^° € £1, there exists a unique tempered solution of (43) with initial 
configuration ^° £ Q.. 

As explained above, the main step to prove this theorem is to obtain an a priori bound 
that we prove in Proposition 3. For a complete proof, we refer to [13] ( or [14, 15]). 

Now we notice that the Gibbs state Vg ^, (j3,A) € (0,+°°) X (— 1,+°°) is formally in- 
variant for the infinite dynamics generated by (§(/))«>o- This can be seen by observing 
that J (Jzf/) (4 )dVp x (d£) = for nice functions / : £2 — >• R. Nevertheless, some care has 
to be taken to prove this. Indeed, we do not know that Jz? is really the generator of the 
semigroup generated by (f)) f >o on the space of bounded measurable functions on £2 in 
the usual Hille-Yosida theory. This can be a very difficult question that we prefer to avoid 
(see [15]). Instead we use the fact that the infinite dynamics can be approximated by finite 
subsystems. 
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Proposition 2. For any j3 > 0,A > — 1, the probability measure Vp x ^ invariant for the 
process (§(f)) t >o. 

Proof. Let n > 2 and consider the local dynamics generated by the generator Jz?„ = srf n + 
where 

^ + (%/)(§) (V« + ^p) 

W)(S) = f -/(I)) 

where / : i2 — > R is a compactly supported continuously differentiable function. The dy- 
namics is essentially finite-dimensional since the particles outside the box {— n — 1 , . . . , n + 
1} are frozen. Thus, the classical Hille-Yosida theory can be applied. The boundary con- 
ditions have been chosen to have 

J(x,fmdv p ^)=o 

for any compactly supported continuously differentiable function / which shows that Vpx 
is invariant for the local dynamics. Since, as a consequence of the a priori bound, the 
infinite dynamics is obtained as a limit of finite local dynamics, this implies that Vp \ is 
invariant for the infinite dynamics. □ 

Then this defines a strongly continuous semigroup of contractions (Pt)t>0 on the Hilbert 
space h 2 (Q.,&,vp i). Moreover, Ito's formula shows that its generator is a closable ex- 
tension of Jzf given by srf + y.y since for any local compactly supported continuously 
differentiable function /, we have 

(p t m)=m+ (\p s xm)ds, ^n, t >o. 

Jo 

A. 1 . Logarithmic energy fluctuation. We have first to consider a clever smooth modifi- 
cation of E. Let < A < 1 and consider a twice continuously differentiable nonincreasing 
function p : R — > (0,1) such that tp(u) = e^ 1 -"' if u> 2, <p(u) = (1 +A + A 2 /2)e~ A if 
u < 1, and <p is concave for u < 3/2, convex if u > 3/2. Finally, < — <p'(«) < X(p(u) < 
e^ l -"\ <p{u) > e~ ;i ( 1 +») and \<p"{u)\ < <p(u) for all u > 0. 
For x £ Z and a > 1 we define the function / as 

f(x,a)= [ cp(\x-y\/a)e- 2 ^dy. 
Jr 

In [15] are proved the following properties on /: 

aexp(-k\x\/o) <f(x,a) < c 2 exp(-A|x|/(7), (44) 

f(x, a) < f{y, o)e 2X \*-y\ , d f(x, a) < e^^daffy, a). (45) 

\d x f(x,a)\ ^minl^^cj),^- 1 /^,^)}, (46) 

g(x)\d x f(x-n,a)\ <4g(\n\ + o)(daf)(x-n,o). (47) 

Here the constants depend only on A . 

For 4 £ (0, +°°) z , fieZ and a > 0, consider the function 

W(§,M ) ff) = £/(*-M.ff)(l + 2&-log&) (48) 
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and let 

W(£) = sup sup a~ l W(£,^,a). (49) 

liez o> g (n) 

Observe that by (44), 

W{Z,n,o)> Cl e- X E(%,n,o), (50) 

for all t, G (0, +°o) z , ju g Z and a > 0. 

For E, e (0, we also consider the function 

^) = sup ^MM . (51) 

The following lemma shows that these two modifications of the logarithmic energy 
fluctuation are equivalent to E. 

Lemma 7. There exists a constant C such that for all £ G (0, +°°) z : 

c-'w(^) < < c#(§), cr 1 ^) < w(^) < c£(t). 

Proof. The inequality W(|) < for all | G (0,+°°) z , is trivial. Let us prove the 

second one by taking a > g(ju), fieZ and denoting 1 + 2% x — log £, x by By (44), we 
have 

W(^fi,a)<c 2 Y J exp(-X\x-fi\/a)H x = c 2 f j e-^ a £ H, 

*EZ n=0 [x-Ml=n 

«=0 |jt-/i|<« 

where the last equality follows from Y,\x~n\=n Hx = L|x-^|<« H x — E|x-/i[<«- 1 ^* and a dis- 
crete integration by parts. Let r > 1 the integer such that r — 1 < g(ju) < r and decompose 
the set {x G Z; |jt — /x| < «} as U^j^Aj where the A ; - are non intersecting intervals of 
length r for j = l,...,K and 1 is of length at most r — 1 . Observe that + 1 is of order 
n/g(jj.). By using (50), we have easily that 

xeAj 

where C depends only on X . Thus we get 

W(£,n, a) < C(l - e- x ' a ) f) e- A "/ CT «W(4) < C'oW{$). 

which concludes the proof of the second inequality. 

The proof of C l E(%) < W(%) < C5(|) for all % G (0,+°o) z , is the same. The first 
inequality follows from (50) and the constant can be taken equal to c\e~^ . The second 
inequality follows from a similar argument to the one used above. □ 

A. 2. The a priori bound. 

Proposition 3 (A priori bound). For each w > 1 there exists a continuous function q w (t), 
t > 0, such that 

p( sup W^(s))>exp(q w (t)g(u)))<e-" 

for each u > 1, t > 0, whenever W(£°) < w and (%{t))t>0 ' s a tempered solution of (43) 
with initial condition 



Proof. We consider a tempered solution (| (f ))r>0 of (43) with initial configuration t, G £2. 
For each £ > 1, jj. G Z and f > we define the stochastic process p^ by 

p t (t)=*B(M)-Cb f g(|Ml + |p*(*)|)Z'(*)^ (52) 



o 



where Q := Co(y, A) is a positive constant that will be chosen later and 

Z(0= fw(^(s))ds. 
Jo 

Since the function /(■) is positive, W'(-) is also positive and this turns Z(-) positive. The 
trajectories of p k are differentiable, decreasing and satisfy P/t+i (?) — Pk{t) < a - s - f° r 
each t > 0. We consider also the sequence of stopping times = inf{f > 0; p k (t) < git-l)} 
which satisfy X k < < +°° and Iirm,^ +O o = +°° a.s. We evaluate now the stochastic 
differential of t — > W(| (f),ju,p/t(f)) for f < (so that p<.(f) > 1). This is given by 

<f[^(f), M ,p*(o)] =i£>{t)dt - c (d a w)(m^,Pk(t))gw+pk(t))w(m)dt 
+ dif\t) 

where 

4 4) (0 = 2 E (/(*-m,p*(0)-/(*+ i-M,P*(0))&(0^-i(0 

' veZ (53) 
+ 52 (/(*+l-M.P*(0)-/(*-l-M.Pt(0))?x(0 

and 

^ = E/^-M,Pfc){2v((^-^_i)^- M ) -v((log^-log^-i)^_ M )}. 

a~gZ 

We first estimate the term Iff' (?) and we show that if Co is taken sufficiently large then, 
for t < Tj we have that 

/«(0 - c (d a w)m^,p k ( t )) g m+p k (t))wm)<o. (54) 

The second term on the right hand side of (53) can be estimated, by using (45) and (46), 
to get to 

f(x+l-fi,p k (t))-f(x-l-ii,Pk(t)) = J da(d x f)(x-fi + a,p k (t)) 

J_da\{d x f){x-p, + a,p k {t))\ 



»i 

< 



< J ^da (d a f)(x-n + a,p k (t)) 
<2 sup d a f(-,p k (t)) 

[x-fi— l^C-/i+l] 

<2e 2X d a f(x-p r ,p k (t)) 



(55) 



which gives us that 

£ (f(x + p k (t)) -f(x-l- /i,P*(t))) Ut) < C £ d(jf{x - m,p*(0)&(0 

A'GZ a-gz 

<C(a a W)(§(0,M,Pt(0)- 
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Now, notice that for any ieZ and for all | G (0,°°) z we have that 

On the other hand, by (50) and since for all x > it holds that log(x) < 1 +x, then we have 
that W(% ,x,g(x)) > c\e^E{<^ ,x,g(x)) > c\e~^% x+ \. Then, we conclude that there exists a 
constant C such that for all x £ Z and % e (0, °°) Z , 

Ui<Cg(x)W(t;). (56) 

To estimate the first term on the right hand side of (53) we use the previous estimate, 
(47) and a similar argument as done in (55). It follows that 

£ (f(x - n,p k (t)) - f(x + 1 - H,p k (t))) Ut)Ui{t) 

xeZ 

< cw{Ut))sM+pk{t)) £ dof(x-nM*))Ut) 

xeZ 

< cw(m)gM+pk(t)) W,m,p*(0)- 

Then, (54) follows. 

(k) 

The term dl\ can be written as 

dl[ k) = £ /(* - P, Pit) {2V ((& - i )<*At_ i ,x) - V ( (log & - log , )rf2V i ,x) } 

agZ 

= £ (/(^ + l-M ) P*)-/(^-P,Pfe)){2V^-Vlo g y^ +1 

xGZ 

= £ (/(x+l-M,P*)-/(x-M,Pfc)){2V^-Vlogy (dN x , x+ i-ydt) 

aGZ 

- 7£ (/(* + 1 - P. Pit) - /(* - P, Pt) ) - V log u A. 

a-gz 

Since the compensated Poisson processes N XJi+ \ (t) — yt are orthogonal martingales with 
quadratic variation y 2 t , then 

<*«? ) = " £ (/(* + 1 -H,Pk)-f{*-H,Pk)) {2V^- Vlog^} (aW„ +1 - y&) 

aGZ 

defines a martingale with a quadratic variation equal to 

J(M«) r = y 2 £ (/(*+ 1 - /i,p t ) -f{x-ii, Pk )f - Vlog£ v } 2 

xGZ 

Using a similar argument to the one in (55), together with the fact that for all x,y € Z 
such that |x|, \y\ < C it holds that \x — y\ 2 < 2C\x — y\, the boundedness of the function /, 
(45), (46), (47) and (56), one has that there exists a constant C such that 

d{Mf ) t <C g m+p k {t))W {£,{!)) d a W{!;{t),pi, Pk {t))dt. 
Similarly we obtain that 



£[f(* + l-M,P*(0)-/(*-M,Pfc(0)]{2V§«W-Vlog&(/)} 
Thus, if the constant Co is chosen sufficiently large, we have 



<Cd a W^(t),H,p k (t)). 



supW(^(f Ax k ),n,p k (t AT k )) < W(%(0),n,kg(n)) +supN(n,k,t) 

t>0 t>0 
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where N{n,k,t)=Mf ) (tAx k )-\ (AfjjPW ■ Observe that exp(Aff } (f A T* ) - ± (Mj® ) ;ATyt ) 
is a martingale with expectation equal to 1 . By the exponential supermartingale inequality, 
we have that 

P(supiV(/i,ik,/) > m) < e~ u . 
t>0 

Thus we proved that for each k > 1, ji g Z and m > 0, 

supW(£(rATt),/i,p t (f AT*)) <W(§(0),/i,*»(M))+« (57) 

r>0 

with a probability greater than 1 — e~ u . Applying (57) for each ft£Z and k > 1 with u 
replaced by u +Akg(p.) where A > 1 is sufficiently large to have £ lt>1 L^gz e~ Akg ^ < 1, 
we obtain 

supW(|(f ATi),/x,p t (f AT*)) < W(§(0),M,*&(M))+^*b(M) + « 
^° (58) 

<kgQi)W(S(p))+Akgbi) + u 

with a probability greater than 1 — e~" uniformly in k and fi. 

Define now k := k t , t > as the smallest integer k > 1 for which p&(/) > g(jli); then 
Tic > t and p^(/) < 2g(fi) as Pk-i(t) < thus choosing = A:, in (58) and using that 

W(<j; ,/x, <t) is increasing in c (since d a f > by the conditions imposed on <p), we get 

"^l'™™ < <^(0)) + A*+-^ <**(§«») -M* +B> 

where in the last inequality we used the fact that g(x) > 1 for all x 6 M. Taking the supre- 
mum over fi and using Lemma 7, we obtain 

f(^(())<cwK(o))+« 

for each t > with probability at least 1 — e _ ". On the other hand, 

2g{pL) > k,g(n) - Co fgQul + \p k (s)\)Z , (s)ds 
Jo 

whence 

k t <2+c [ tg M+\pm Z ' (s)ds . 

Jo g(p) 

Since p k (s) < k t g(jj.) for any s S [0,f] and g is increasing, we have that g(\ji\ + \Pk(s)\) < 
g(ju + ktg(n)). On the other hand for x > 2, g(x) < x together with the fact that for x,y G 
^ g(Wbl) < <?(bl)<?(bl) an d since g(l + x) < 1 + g(x) for x > 1, we obtain that g{jx + 
kig(jj,)) < g(p.)(l + g(k t )). As a consequence we obtain that 

k t <2 + C Q Z(t)(l+g(k,)). (59) 

Since for all x > 1 we have that g(x) < 1 +2 y^jxf, then 



^,<2 + C Z(f)(2 + 2Vfe)- 
Finally, it follows that < 2 + 4CoZ(f ). Then, since g is increasing and by plugging the 
previous inequality in (59), we obtain that 

k, < 2 + C Z(f)(l +g((2 + 4CoZ{t)) 2 )). 

Recalling that Z'(f ) = (f )) we obtained that there exists a constant M > depending 
only on A such that for any w > 1 and any initial condition ^ (0) satisfying (0)) < w, 



supM~ l Z'(t)-w{l+Z(t)g(Z(t))) < u 



> 1 -< 
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The a priori bound follows from this last inequality (see [13], Proposition 1). □ 
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